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The X-ray afterglow of GRB030329 
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Abstract. We report on XMM-Newton and Rossi-XTE observations of the bright (fluence ~ 10 erg cm ) and 
nearby (z=0.1685) Gamma- Ray Burst GRB030329 associated to SN2003dh. The first Rossi-XTE observation, 5 
hours after the burst, shows a flux decreasing with time as a power law with index 0.9±0.3. Such a decay law is 
only marginally consistent with a further Rossi-XTE measurement (at t-toRB ~30 hr). Late time observations 
of this bright afterglow at X-ray wavelengths have the advantage, compared to optical observations, of not being 
affected by contributions from the supernova and host galaxy. A first XMM-Newton observation, at t-tGRB ~37 
days, shows a flux of 4xl0 -14 erg cm -2 s _1 (0.2-10 keV). The spectrum is a power law with photon index 
r=1.9 and absorption <2.5xl0 20 cm -2 , consistent with the Galactic value. A further XMM-Newton pointing at 
t-tGRB ~61 days shows a flux fainter by a factor ~2. The combined Rossi-XTE and XMM-Newton measurements 
require a break at t~0.5 days in the afterglow decay, with a power law index increasing from 0.9 to 1.9, similar 
to what is observed in the early part of the optical afterglow. The extrapolation of the XMM-Newton spectra to 
optical frequencies lies a factor of ~ 10 below simultaneous measurements. This is likely due to the presence of 
SN2003dh. 
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1. Introduction 

A very bright Gamma-Ray Burst (GRB) has been 
recorded by several satellites on March 29, 2003. The 
accurate localization obtained with HETE-2 after about 
1 hour jVanderspek et al. 2003| IRicker 20 03) prompted 
rapid observations that revealed a bright optical tran- 
sient with R magnitude about 13 ({Peterson fc Price~2 003 
ITorii 2003|) . A redshift of z=0.1685 has been measured 
for the GRB030329 host galaxy ljGreiner et al. 20031 
IGaldwell et al. 2003}) . This is the second smallest red- 
shift determined for a GRB (the smallest one is that of 
GRB980425, likely as sociated with SN1998bw at z=0.0085 
IjGalama et al. 1998)) as supported by recent observations 
(Pia n et al. 20 03 )). Although the intrinsic luminosity of 
GRB030329 was in the low end of the distribution for 
GRBs (see below), its proximity led to a very high flu- 
ence for the prompt emission (1.2xl0 -4 erg cm -2 , 30-400 
keV, IRicker 2003fl and to a bright afterglow which can 
be studied with unprecedented detail at all wavelengths 
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and over long timescales after the burst explosion. Among 
the most interesting results reported so far, is the finding 
of clear spectroscopic signatures of an underlying super- 
nova (Stanek et al. (2003a), Hjorth et al. (2003)), which 
provides strong evidence for the association between long 
GRBs and core collapse supernovae. 

The early phases of the X-ray afterglow of 
GRB030329 were observed with two Rossi-XTE point- 
ings obtained 5 hours and 1.24 days after the 
burst ({Marshall fc Swank" 120031 Marshall et al. 2003|l . 
Unfortunately, no further X-ray data could be collected 
during the following month. At the beginning of May 
the GRB position became compatible with the visibility 
constraints of the XMM-Newton satellite, which per- 
formed two observations 37 days and 61 days after the 
burst. The brightness of GRB030329 , coupled with the 
large collecting area of XMM-Newton , allow us to study 
in detail for the first time an X-ray afterglow at such 
long times after the prompt emission. This is particularly 
interesting since at these wavelengths we do not expect 
significant contamination from the underlying supernova 
and/or host galaxy. 
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2. Data analysis and results 

2.1. RossiXTE 

The first Rossi-XTE observation, consisting of two time 
intervals of duration 1500 s and 500 s, respectively, was 
done on March 29, about 5 hours after the GRB. The 
Proportional Counter Array instrument (PCA, Jahoda et 
al. 1996) aboard Rossi-XTE consists of five Proportional 
Counter Units (PCU). Only three (n. 0,2 and 3) and two 
(n.O and 2) of them were on during the first and second 
time interval, respectively. Since the instrumental gain is 
known to vary between PCUs, we separately extracted the 
spectra from the two time intervals. In order to increase 
the signal to noise ratio, only the top layer anodes were 
used in the analysis. The corresponding response matrices 
were generated using PCARSP V8.0 and the background 
spectra were estimated with the faint-source model as in- 
put to PCABACKEST V3.0. All the errors quoted below 
are at the 90% confidence level. 

Since the source spectral shape did not vary between 
the two intervals (except for the normalization), we fit- 
ted them together, obtaining a best fit with a power 
law model with photon index r=2.17ioQ3 an d absorption 
Njj <5xl0 21 cm~ 2 . The average flux in the first inter- 
val was F x =(1.38±^)xl0- 10 erg cuT 2 s" 1 (2-10 keV). 
During this observation the flux decreased with time as a 
power law with index <5=0.9±0.3 (see Fig.0. 

A second observation was performed 30 hours after 
the GRB, on March 30 from 17:30 to 19 UT, with a net 
observing time of 3500 s. Due to the lower flux, the derived 
spectral parameters are less constrained. A power law fit 
yields r=1.8±g;l, N H <3.7xl0 22 cm- 2 , and F x =(lt°; 2 ) x 
lO" 11 erg cm" 2 s" 1 (2-10 keV). This flux is smaller than 
the extrapolation of the power law decay found in the first 
observation, suggesting the presence of a break in the light 
curve. A single power law fit to all the Rossi-XTE data 
gives only a marginally acceptable fit for a slope 8 ~1.5 
(x 2 =7.6 for 3 degrees of freedom (dof)). 

Finally, we analyzed an observation carried out on 
April 6. The afterglow was not detected, with a flux upper 
limit of 2.6 xl0~ 12 erg cm" 2 s _1 (3er). 

2.2. XMM-Newton 

XMM-Newton observed the position of GRB030329 start- 
ing on May 5, 2003 at 12:30 UT, for an observation 
length of ^12 hours. The last part of the observation 
was affected by high particle background and was there- 
fore excluded from our analysis, resulting in net exposure 
times of 29 and 32 ks, respectively in the PN and MOS 
cameras of the EPIC instrument iStriider et al. 2 001 
ITurner et al. 200 All the cameras operated in Full 
Frame mode and with the thin optical blocking filter. The 
data were processed using SAS version 5.4.1. 

A source with a PN net count rate of 0.018 counts 
s" 1 was detected at R.A. = 10 h 44 m 49.9 s , Dec. = +21° 
31' 15" (J2000, error radius of 4"), consistent with the 
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Fig. 1. X-ray afterglow of GRB030329 as measured dur- 
ing the first Rossi-XTE observation. Each time bin is 500 
s long. The line is the best fit with a power law of index 
<5=0.9±0.3. 

position of GRB030329 . Its flux during the observation 
is consistent with a constant value (see Fig- - X second 
source of similar intensity (0.022 PN counts s _1 ) is present 
at a distance of ~30" to the NW. Its spectrum (a power 
law with photon index T=1.7±0.2 and N H =(5±2) x 10 21 
cm~ 2 ) and its positional coincidence with a galaxy at 
z=0.136 (Kri sciunas et al. 2003J1 indicate that this source 
is an AGN. 

To measure the afterglow spectrum we used an extrac- 
tion radius of 15" in order to minimize the contamination 
from the AGN. For the extraction of the background spec- 
trum we chose a circular region (radius 15") at the same 
distance from the AGN as the afterglow. This was done 
to take into account the small contamination from the 
AGN to the spectrum (we estimate that at most 20% of 
the counts in the source extraction region could be due 
to the AGN). The spectra, over the 0.2-10 keV energy 
range, were rebinned in order to have at least 30 counts 
per channel. After checking that consistent results were 
obtained in the three cameras, we fitted jointly the MOS 
and PN data. 

The best fit (x 2 /do/=16.5/15, see Fig. EJ) was ob- 
tained with an absorbed power law with photon in- 
dex T = 1.92t°- 2 ^, N H <2.5xl0 20 cm" 2 , and flux 
F x =(4.2±0.5)xl0~ 14 erg cm- 2 s" 1 (0.2-10 keV). The 
Galactic absorption in this direction is Nh=2x10 20 
cm~ 2 (Dick ey~fc Lockman 1990] ). Other simple models, 
e.g. blackbody, thermal bremsstrahlung, thermal plasma 
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Fig. 2. Background subtracted PN light curve of the X- 
ray afterglow during the first (left) and second (right) 
XMM-Newton observation. The bin sizes are 2,000s and 
10,000 s, respectively. Note that the count rate quoted in 
the text is higher than the one shown here since it has been 
corrected for the fraction of source counts falling outside 
the extraction region. 



(MEKAL) gave unacceptable fits (the respective values of 
X 2 /dof being 73.5/15, 26.3/15 and 36.8/15). 




Energy (keV) 

Fig. 3. EPIC spectrum of the X-ray afterglow of 
GRB030329 fitted with a power law model. Upper line 
and data refer to the PN, lower ones to the MOS. The 
bottom panel shows the best fit residuals in units of stan- 
dard deviations. 



A second XMM-Newton pointing started on May 28, 
2003 at 21:00 UT. It lasted about one day, but it was 
severely affected by periods of high particle background, 
which were excluded in our analysis, resulting in about 
40 ks of useful data. The analysis was performed as de- 
scribed above for the first observation. The source at the 
GRB030329 position fainted to F x =(2.2±° j) x 10~ 14 erg 
cm~ 2 s _1 (0.2-10 keV), thus confirming that it is indeed 
the GRB afterglow. Its spectrum could be fit by an ab- 
sorbed power law with T = 2.lt°[2 and n h <4.5x10 20 
cm~ 2 (x 2 /dof=4.5/13). Also in this case a thermal model 
gave a worse fit (x 2 /dof=21.3/13 for a MEKAL). 

3. Discussion 

All our measurements of the X-ray afterglow of 
GRB030329 for the 2-10 keV range are plotted in Fig. 

where one can see that the break in the afterglow decay 
suggested by the second Rossi-XTE observation is clearly 
confirmed by XMM-Newton . While a power law index 
5=0.9 was found during the first Rossi-XTE observation, 
a fit to the following points gives a slope <5=1.86±0.06. We 
estimate that the break occurred in the time interval 0.3 
- 0.8 days, with a most likely value of tbreak=0.45 days. 
This value is consistent with the time of the break in the 
optical afterglow (tbroak=0.48 days, Price et al. 2003). 

This achromatic break can be readily explained as a 
"jet-break", due to the decreasing bulk Lorentz factor T, 
making 1/r equal to the jet opening angle (see Rhoads 
1999). Following Frail et al. (2001) it is then possible to 
estimate the opening angle of the jet (~ 3°), correspond- 
ing to r ~ 19 at the time of the break. With this opening 
angle, the "true" energy radiated by the burst in 7-rays 
turns out to be E y = 3 x 10 49 erg (see Frail et al. 2001 
for the relevant uncertainties concerning these estimates) . 
This value is at the very low end of the distribution found 
by Frail et al. (2001), making GRB030329 an atypically 
weak burst. Note also that the optical light curve shows 
several achromatic "rebrightnenings" and breaks (Granot 
et al. 2003, and references therein), which weaken the asso- 
ciation of the first break with the jet-break. The paucity 
of the X-ray data does not allow us to infer if the X- 
ray follows the optical during the several rebrightenings 
occurring in the optical band, which could help to investi- 
gate the origin of such rebrightenings (Lazzati et al. 2002; 
Granot et al. 2003). 

In Fig. [S] we show the simultaneous optical and X-ray 
spectra corresponding to the second Rossi-XTE and to 
the first XMM-Newton observation epochs. Since we did 
not find an exactly simultaneous spectrum at ~ 30 hours 
after the trigger, we show in Fig.[S]the optical points taken 
immediately before (22 hours) and after (40 hours) the 
Rossi-XTE observations. As can be seen, at this epoch the 
extrapolation of the optical spectrum joins very smoothly 
the X-ray data, and fits both their normalization and 
slope. Optical and X-ray fluxes therefore belong to the 
same spectral segment characterized by F(v) oc v^ 1 (i.e. 
a flat spectrum in vF(v)). This shape can be explained by 
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the standard synchrotron-external shock model (e.g. Sari 
et al. 1998) as due to a population of relativistic electrons 
injected in the emitting region with an energy distribution 
oc 7~ p with p ~ 2, with both the optical and X-ray fre- 
quencies laying (at ^30 hr) above the cooling frequency 
v c (this is the frequency produced by those electrons that 
have just cooled in a dynamical time). 

The optical and X-ray decay slopes before and imme- 
diately after the break at ~ 0.5 days are the same (i.e. 
F(t) oc t~ - 9 and F(t) oc i -1 - 9 ). The first decay slope is 
consistent with what expected in the case of p ~ 2, slow 
cooling regime and the cooling frequency below the opti- 
cal (Panaitescu & Kumar, 2000). In this case the decay 
index is independent on the density profile. The second 
decay slope is instead consistent with what expected after 
the jet break (t~ p ) if the jet matter expands laterally at 
a velocity close to the speed of light (Rhoads 1999). 

The evolution of the high-energy spectrum depends 
on the behavior in time of the cooling frequency which 
in turn depends on the circumburst density profile and 
on the jet dynamics. Therefore, after 30 hours, there are 
two possibilities. If the circumburst medium is homoge- 
neous, the cooling frequency v c decreases as i -1 / 2 (before 
the jet break) or remains constant (after the jet break with 
sideway expansion) , leaving the optical to X-ray slope un- 
altered. Instead, in the case of a R~ 2 wind density profile, 
the corresponding behaviours are v c oc t 1 ^ 2 and i/ c =const. 
Only in the wind case before the jet break v c possibly 
overtakes the optical band (see e.g. Panaitescu & Kumar 
2000). This gives a spectral break between the optical and 
the X-ray bands, accompanied by a flattening of the opti- 
cal light curve (due to the appearance of the v < v c harder 
spectral slope in the optical band). We conclude that the 
optical to X-ray flux ratio of the synchrotron-external 
shock component either remains fixed or decreases in time. 

The optical data of May 5th define a steep (spectral 
index a opt > 1) spectrum and lay above the extrapolation 
of the X-ray spectrum. We interpret this as evidence for a 
significant contribution, in the optical, of SN2003dh. The 
X-ray flux can be used to estimate an upper limit to the 
flux of the optical synchrotron-external shock component 
(indicated by the dotted line in Fig. EJ: vF(v) < 10~ 14 
erg cm~ 2 s _1 . It is an upper limit if v c > v pti in the case 
of homogeneous circumburst density this becomes the ac- 
tual value. We conclude that the optical lightcurve, around 
May 5th, is dominated by the light from the supernova 
SN2003dh, by at least two magnitudes (see also Fig. 0}. 

Similar considerations apply to the data at t-tcRB = 
61 days, although with larger uncertainties owing to the 
less constrained XMM-Newton spectral slope. 

4. Summary 

Thanks to the high sensitivity of XMM-Newton we 
could study the optical-X-ray SED of the afterglow of 
GRB030329 and its time evolution up to late times. This 
is particularly important for this burst due to its associa- 
tion with the supernova 2003dh, which, at late times, con- 
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Fig. 4. The X-ray light curve of GRB030329 (large 
dots and upper limit at ~8 days) is compared with 
the optical data at times close to the epochs of the 
X-ray observations. Optical data are from Burenin et 
al. (2003a); Burenin et al. (2003b); Fitzgerald & Orosz 
(2003); Ibrahimov et al. (2003); Price & Mattei (2003); 
Price (2003); Rykoff & Smith (2003); Stanek et al. 
(2003b); Stanek et al. (2003c); Zharikov et al. (2003). 
The upper limit at 37 days corresponds to the optical 
flux calculated through the extrapolation of the XMM- 
Newton spetrum (see Fig. [SJ . The dotted lines are only 
indicative of the time decay slopes, and are not fits to the 
optical data. 

tributes to the optical flux (Stanek et al. (2003a), Hjorth 
et al. (2003)). The early and late time X-ray data, com- 
bined with the simultaneous optical detections, have been 
used to estimate the contributions of the (non-thermal) 
afterglow and supernova components at optical frequen- 
cies. 

Our main results are the following: 

— The first epoch Rossi-XTE data define an X-ray light 
curve decaying in time as t~ 9 . This decay index is 
consistent with the one of the optical flux. 

— The two XMM-Newton observations at late epochs 
yield spectra well fit by a power law with photon index 
~2, 

— The Rossi-XTE and XMM-Newton data, taken to- 
gether, are consistent with a break in the light curve 
occurring at ~0.5 days, simultaneously with the op- 
tical break. After this break the afterglow decays as 
t- 1 - 9 . 

— The optical to X-ray SED at 30 hours strongly indi- 
cates that both spectral bands lay on the same branch, 
above the cooling frequency. 
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Fig. 5. The quasi simultaneous optical to X-ray SED of 
GRB030329 at the epochs of the second Rossi-XTE and 
of the first XMM-Newton observation. Optical data are 
from Zharikov et al. (2003), Fitzgerald & Orosz (2003) and 
Ibrahimov et al. (2003). Open circles refer to fluxes de- 
reddened assuming an extinction of Ay = 0.16 consistent 
with both the Galactic value of the column density and 
the Nh found for the XMM-Newton fit. Triangles assume 
instead Ay — 0. For the SED at ~ 30 hours, the optical 
and X-ray data lays on the same power law [F(v) oc v , 
top dotted line]. The bottom dotted line corresponds to 
the same spectral shape for the May 5th SED. 

- The optical to X-ray SED on May 5*' 1 (and possibly 
also on May 28" 1 ) indicates instead an optical excess 
that we interpret as due to SN 2003dh, which should 
dominate (by a factor ^10) the non-thermal optical 
emission. 
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